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Abstract 



We give a selection of results on spectrum and decay constants of light and 
heavy-light hadrons. Effective fields theories relevant for their lattice calculation, 
namely non-relativistic QCD (NRQCD) for heavy quarks on the lattice and Chiral 
I Perturbation Theory for light quarks, are briefly discussed. 

^ 1 INTRODUCTION 
o 

O I The Standard Model of the strong and electroweak interactions is based on a SU{3) x 
SU{2) X U{1) gauge symmetry with three generations of quarks and leptons as fermionic 
matter fields and a scalar field, the Higgs, which is responsible for the masses of the weak 
Q ■ SU{2) gauge bosons and of the fermions. For a recent review about the status of the 
Standard Model and new physics see e.g. [1]. 
■ The SU{3) 'sector' of the model is Quantum Chromodynamics (QCD), a gauge theory 

of the strong interaction. With relativistic Dirac quarks, the model can be described 
classically by the Lagrangian 
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Cqcd = q{i^,D^-m)q--F;,F^^''. (1) 



^ \ The q fields are 4-component Dirac spinors, and the are covariant derivatives, e.g. 

= — igsA^j^f^ with [D^, D^] = —iggF^j^t^, where are the field strength ten- 
sors, Qs is the coupling constant, and the are generators of SU{?>) in the fundamental 
representation. A consequence of the self-interactions among the gluon fields is asymp- 
totic freedom, i.e. the interactions between particles become weak at short distances and 
can be described with perturbation theory in the strong coupling as = g1/{4:Ti). At larger 
distance, the forces become strong, and non-perturbative methods are necessary to under- 
stand how hadron masses arise and whether it is possible to explain the hadron spectrum 
from first principles within the theory of strong interactions. 

Although rather successful, the Standard Model by itself does not seem completely 
satisfactory. On the experimental side, recent discoveries such as neutrino mixings, new 
results from accelerator experiments [3,4] and indications for 'dark energy' in the cosmos 
indicate a need for an extension of the model. The Higgs particle has not yet been 
found; recent reviews of the status of Higgs searches are [5] . Further there are theoretical 
motivations to search for physics beyond the Standard Model (for a discussion see e.g. [2]). 
The Standard Model contains a considerably large set of coupling constants and masses 
as input parameters. It does not explain the values of typical energy scales such as the 
masses of the weak gauge bosons. 

A strategy in the research is to simultaneously measure as many physical quantities 
as possible, test the results for self-consistency within the Standard Model and search for 
indications of new physics. Among the most interesting search grounds are the elements of 
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the Cabibbo-Kobayashi-Maskawa (CKM) matrix which parameterizes the flavor changing 
weak currents and provides a mechanism for CP violation within the Standard Model. 
Those CKM matrix elements which are relevant to reactions of heavy, for example b and 
c, quarks are at present studied intensively in experiment and theory. We introduce the 
CKM matrix with an emphasis on B meson decays in the framework of the weak effective 
theory in Section 11.11 The status of the CKM matrix is reviewed in [6] . For a review 
about recent results on quark masses see Ref. [7]. 

Description of the long-range interactions of QCD requires non-perturbative tech- 
niques. Using a four-dimensional lattice description of space and time it is possible to 
calculate matrix elements numerically on a computer within a path integral formalism. A 
brief introduction to the lattice formalism is given in Section |2j for detailed recent reviews 
see [8]. 

Ideally, the lattice extent L should be much larger than the extent or the Compton 
wavelength of the particles that are supposed to be described, and the inverse lattice 
spacing a should be much larger than the masses and momenta in the theory in order to 
avoid cutoff effects. The lightest hadrons, the pions, have a mass of around 140 MeV, 
whereas the B meson has a mass of 5.28 GeV and contains a heavy quark with a mass 
of 5 GeV. The problem is how lattice simulations can accommodate this large range of 
scales. 

Ideally <C masses and <C 

energy splittings 

In Reality 

L = 2 - 3 fm = 2-4 GeV 

(L-i = 0.07-0.1 GeV) 

finite size effects cutoff effects 

To calculate properties of hadrons with b quarks on the lattice, one can for example 
simulate at lighter quark masses where discretization errors are under better control, and 
use extrapolations in the heavy quark mass. Fortunately the energy level splittings of 
b hadrons are much smaller than their masses: of the order of A.qcd = 200 — 500 MeV 
or smaller, where Aqcd is the energy scale where QCD becomes non-perturbative. The 
dynamics of the heavy quarks can be accounted for as small corrections proportional to 
powers of the inverse heavy quark mass. This is the basis for effective field theories devel- 
oped for heavy quarks: Heavy Quark Effective Theory (HQET) [9] and non-relativistic 
QCD (NRQCD) [11,12] (for reviews see [10]), which can be used to simulate heavy quarks 
directly on the lattice while avoiding large discretization errors due to the large mass. We 
discuss NRQCD in Sectional 

Practical simulations with light quarks are computationally expensive and sensitive to 
the finite lattice volume. Therefore one often uses quark masses much heavier than u and 
d quark masses and extrapolates the results to the physical values of the quark masses. A 
formalism for this can be derived using chiral perturbation theory (xPT), an expansion 
around the chiral (zero quark mass) limit describing low-ener gy d egrees of freedom of 
QCD such as pions and nucleons. This is introduced in Section 

Lattice results for light and heavy-light hadron masses and heavy-light current matrix 
elements are discussed in Section EJ 



2 



1.1 Heavy quark decays 



Study of weak decays of quarks is of interest for determinations of elements of the CKM 
matrix which parameterizes the mixings of quark generations in the Standard Model: 

£ = (uL-CLtL) nv,,,,) (^t^ Wl + h.c. , (2) 

where ul, di, cl, sl, ti, are left-handed quark spinors, VF^ a charged weak gauge boson 
and (?2 the weak gauge coupling. V is a unitary matrix. There are indications that some 
of its elements have a non-trivial complex phase giving rise to CP violation. The CKM 
matrix elements with presently the largest uncertainties are the ones relevant to decays 
or mixings of the b quark: Vet, Vub and Vtd- iKifeP describes for example the leptonic 
meson decay I'^J^i, where / is a lepton (e, /.i or r) and ui the corresponding neutrino 

and semileptonic decays into a light meson {B vr, p, u) and a lepton-neutrino pair, 
and iKfeP determines the semileptonic decay B ^ D and a lepton-neutrino pair. |V^dp 

is proportional to the oscillation frequency between the mass eig enstates of the 5° - B^ 
mixing, which is described by the left and middle diagrams in Figure Q in the electroweak 
theory. Processes at energy scales much less than the W boson mass can be calculated 




Figure 1: Box diagrams describing B^ — i?" mixing in the electroweak theory (left and 
middle) and the weak effective theory (right). 

within the weak effective theory where interactions mediated by the W or Z particles can 
be described by point interactions. B^ — B^ is described by the third diagram in Fig. 

To relate the weak processes between quarks with exclusive reaction rates of mesons, 
one uses form factors which get contributions from long-distance QCD interactions, and 
therefore have to be calculated nonperturbatively. This can be done from first principles 
using the lattice. In the effective theory the B meson decay is described by a matrix 
element of the heavy-light axial vector current 

(0|A^(x)|5(p))=^/BP;.e-*^^ (3) 
where Jb is the B decay constant. The branching ratio for the decay B~^ l^vi is 

BR(B-^l-.,) = %£!!i(l-4)V||y„,|V (4) 

Stt V i^bJ 

where Gp = g2/{^^w) the Fermi constant and tb the B lifetime. If fs is known, 
\Vub\ can in principle be determined experimentally from this decay. There exist only 
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experimental upper bounds on fs, fss ^iid /d, but there are results on fn/- 

fn, = 280(17) (25) (35) MeV [13] and 
= 285(19)(40) MeV [14]. 

In the weak effective theory, the form factor for the — mixing matrix element can 
be parameterized as fsBs, where the "bag parameter" Bb quantifies to what extent the 
matrix element is described by 5-to-vacuum currents: 

^ _ 3 mihi,dL)ihi,dL)\B<^) 

2{B^\h^^,d\0){0\ln,l,d\B^)' 

The oscillation frequency of the mass eigenstates is proportional to the mass difference 
and related to the form factors by 

AMd (X \V,lVtd\^flBB. (6) 



2 QCD ON THE LATTICE 



2.1 Gauge fields 

Matter fields, e.g. quarks, sit on the lattice sites which are separated by a spacing a. 
The gauge fields on the lattice are represented by fields U^,^ = exp[igsa f^_^^dzAlj^{z)t'^] = 
exp[—igsaA1j^{x)t'^], parallel transporters between neighboring lattice sites. Line integrals 
of gauge fields over closed paths are called Wilson loops. The smallest (1x1) Wilson loop 
on the lattice is a product of gauge links over the the nearest neighbors, called plaquette. 



Ux,fiUx+^i,uU^^jjJJl. 



By expanding around a ^ 0, one finds a tree- level relation between the plaquette and 
the continuum field strengths. Thereby one obtains the Wilson (or 'plaquette') lattice 
action as a discretization of the continuum gauge field action of Eq. (P), 



□ 



2iV, 



Ny. number of colors, 



(7) 



which has lattice spacing errors at O(a^). To further reduce discretization effects, actions 
can be improved. For gauge field actions, this consists of adding larger Wilson loops 
(1 X 2, ...). Improvement can be done by removing discretization effects order by order 
in a (originally suggested by Symanzik for scalar field theory [15], and developed into an 
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improvement program for on-shell quantities in QCD by Ref. [16]), or with renormalization 
group methods to obtain renormahzation group improved (RG) 17] or perfect [18] actions. 

At typical values of (3 in lattice simulations, there are large corrections due to gauge 
field loops on the lattice which shift the expectation value of U substantially with respect 
to the free field value, one. The perturbative corrections can be reduced with 'mean- 
field' (or 'tadpole') improvement [12]: the gauge links U are divided by their expectation 
value which can be calculated in perturbation theory or determined nonperturbatively in 
simulations. 

2.2 Lattice fermions 

Discretization of the Euclidean continuum Dirac action by substituting the covariant 
derivatives by covariant symmetric lattice differences gives the 'naive' lattice fermion 
action 



which is chirally symmetric if m — * and has O(a^) errors. However, the naive discretiza- 
tion leads to a flavor multiplication, the so-called 'doublers'. At m = 0, the free fermion 
propagator has a pole at fc^ = as in the continuum but also poles at fc^ = vr/a. There 
are 16 species of fermions, which occur in pairs of opposite chirality. 

Wilson's solution to the doubling problem is to add a term of the form J2x 'Mx^lx to 
the action, where A is a covariant second derivative. The doublers obtain masses which 
remain finite in lattice units: ma 7^ if a 0. Chiral symmetry receives corrections 
at a 7^ 0, and 0{am),0{ap) discretization errors occur. 0(a) errors can be removed 
from the action with the clover term proportional to a'^'^x'^x'^iJ.i^G^'^Qx, where G^u is 
a discretized version of the field strength tensor using four neighboring plaquettes [19]. 
The coefficient of the clover term can be calculated in perturbation theory (a common 
choice is at tree- level using tadpole- improvement). Most recent calculations use a non- 
perturbative determination of the clover coefficient [20] and are 0{a) improved to all 
orders in perturbation theory. 

Staggered fermions are obtained from a spin-diagonalization of naive fermions: 



with 7^. = 7i^72^73^74*, where x = {xi, X2, X3, X4). The x fields are one- component. In the 
massless case, the theory has a U{1) x U{1) chiral symmetry at finite a and a U{A) x U{A) 
chiral symmetry in the continuum limit. Discretization errors are O(a^). Improvement is 
possible by adding higher dimensional operators. 

Lattice formalisms for doubler-free, chiral fermions are given in [21-23]. 

2.3 Extracting physical quantities 

Green functions can be calculated by evaluating the path integral over the lattice degrees 
of freedom numerically. For example, a two-point function of a field 0{x) which can be 
composed of more elementary fields {0j(x)} is given by 





Qx IxXx ; Qx Xxlx 5 



(9) 




(10) 
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where Vcf) denotes integration over all dynamical fields (gauge, fermion, etc..) in the 
theory. To determine for example fs from the lattice, it is necessary to calculate the 
renormalization factors to match the unrenormalized lattice matrix element of the axial 
vector current to the corresponding matrix element in continuum QCD. 

Ideally, these calculations are done at various values of the lattice spacings, and the 
continuum estimate is obtained by extrapolating as a function of a to a — 0. In practice, 
some lattice calculations are performed only at one or two values of a, in which case 
a continuum limit cannot be taken, and the discretization effects have to be included 
into the estimate of systematic errors. With NRQ CD calculations, higher dimensional 
operators are included as discussed in Section Kllf and an a — >■ extrapolation cannot 
be done out of principle. Calculation at several values of a then serves to determine the 
systematic error from keeping the lattice spacing finite. 

In full QCD, the path integral includes gauge and fermionic fields j VU'Dq'Dq. To 
decrease computational expenses, many calculations are done in the quenched approxi- 
mation, i.e. the vacuum polarization due to quark loops is neglected. 

To use lattice results in phenomenology, it is necessary to estimate systematical errors 
as accurately as possible. The most important sources are: 

• Finite lattice spacing 

• Finiteness of lattice volume 

• Quenching (unphysical number of dynamical quarks) 

3 EFFECTIVE THEORIES AND THE LATTICE 
3.1 NRQCD 

Non-relativistic QCD (NRQCD) is an effective theory formulated for heavy quarks as- 
suming that their dynamics is non-relativistic, with correction terms which can be added 
within a systematic expansion. For quarkonia the higher order interactions are arranged 
in a expansion, where v is the heavy quark velocity (see e.g. [12]). In heavy-light sys- 
tems it is an expansion in v or 1/M, where M is the heavy quark mass. At infinite mass, 
the heavy quark is just a source of the color electromagnetic field, whereas at finite M, 
there is a recoil of the heavy quark due to the interaction with soft gluons with typical 
momenta of 0{Kqcd)- One can argue that the heavy quark momentum Pq and light 
quark momentum pq are equal due to momentum conservation within the rest frame of 
the meson: 

Mv-Pq = Pq^OiAgcD). (11) 

Therefore v ~ Aqcd/M ~ 0.1 in S mesons and should be a reasonable expansion param- 
eter to specify corrections to the M — s> oo (static) limit. Contributions at 0{1/M) are 
the kinetic and the spin-colormagnetic energy of the heavy quark; at 0(1/M^) a heavy 
quark spin-orbit interaction and a Darwin term are added. The 0(1/M^) Lagrangian for 
the heavy quark is given by 

C = ^\Dt + H)i;, (12) 
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with heavy quark Pauh spinor ifj and the Hamiltonian 



a-B 



H 



2M 



+ 



igsC2 

8M2 



(D-E-E-D) 



gsCs 
8M2 



a - {D X E-E X D) 




(13) 



8M3 



The last term is the first relativistic correction to the kinetic energy of the heavy quark, 
which is usually included in calculations at 0(1/M^). The coefficients of the various terms 
can be found with matching calculations to full QCD in the continuum. In the lattice 
calculations described in Section HJ they are set to their tree-level value one, using mean 
field improved gauge links. 

Discretizing the Lagrangian (fT^ one can simulate b quarks directly on the lattice, 
since there are no 0((aM)") discretization errors. Errors 0{o?fP') and 0{aMv'^ /2) arising 
from discretization of the spatial and temporal derivatives in the NRQCD Lagrangian 
can be corrected for by adding further terms to the Hamiltonian (fT!^ . If matrix elements 
of operators are to be calculated in this formalism, the 1/M corrections to the operators 
have to be taken into account as well. Simulation of the 1/M corrections within HQET 
on the lattice using nonperturbative renormalization is discussed in [24]. 

Other methods to avoid large discretization errors used in the calculations discussed 
here are to simulate heavy quarks around the charm and extrapolate to the b, or to use a 
non-relativistic interpretation of a Wilson or clover heavy quark action called FNAL [25] 
in this article. Refs. [25] and [26] formulate on-shell improvement programs for heavy 
quarks by adding further operators to the Wilson or clover action. For recent reviews see 
also [27,28]. 

3.2 Chiral perturbation theory 

The chiral symmetry of massless QCD can be understood as spontaneously broken. As 
a consequence one would expect massless Goldstone bosons in the spectrum. One can 
identify the physical pions with the Goldstone bosons if small quark mass terms are 
included in the QCD Lagrangian resulting in a small explicit breaking of chiral symmetry 
and a finite but small pion mass m^. Typical momenta of low-energy interactions of 
pions will be O(m^). Chiral perturbation theory can be formulated as an effective theory 
expanding in powers of the pion mass and external momenta p = 0{mT^), the so-called 
p expansion [29]. It is used to describe physical interactions at low energy scales, for 
example pion-nucleon scattering. It gives predictions for the expansion of hadron masses 
around the zero quark mass limit, which can be used in the analysis of lattice calculations 
to extrapolate lattice hadron masses simulated at larger quark mass to the physical light 
{u,d) quark mass. The leading dependence at O(p^) is expected to be linear in the 
quark mass rriq where rrig oc m^, with corrections at 0{p^) etc. To facilitate inclusion 
of interactions with baryons within the p expansion. Heavy Baryon Chiral Perturbation 
Theory (HBxPT) was developed [30], which works in the limit of infinite baryon mass. 

The quark masses in present simulations are such that pions are not much lighter than 
the lattice nucleons, and a formulation with relativistic nucleons seems more appropri- 
ate [31,32]. 

For the pions one can use a representation based on their nature as Goldstone bosons 
of the broken SU{2) x SU (2) chiral symmetry of massless QCD with two flavors of u and 
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d quarks with 

= -Jl Tr [d.U^d^U + x^U + xU^)] ■ (14) 

U = exp(-^T7f) transforms according to the (2,2) representation of SU(2) x SU(2). x = 

2BoAi, Ai is the quark mass matrix and Bq is proportional to the chiral condensate. The 
7? are the pion fields, and the Tj are Pauli matrices, /^r is the pion decay constant. 
The nucleon Lagrangian at lowest order (0(p^)) is 

CW = ^ i^t^^D'^ - mo) * + ^(?^^7/.75W^*, (15) 

with V? = U, Dfj, = + 1[m"'',9^m] and = iu^d^Uu^. \E' is the Dirac spinor of the 
nucleon, mo the nucleon mass in the chiral limit, and Qa the nucleon axial vector coupling 
in the chiral limit. The 0{p'^) Lagrangian is given by 

£(2) = ciTr(x+)^^ + -^TiiuaU^ f^D^D^^I/ + h.c.) 

+ |Tr(MX)^*- J^7^71wm'«-]*> (16) 

with x+ = M^xM^ + ux^u. Tr() refers to the trace over the flavor indices. 

Chiral extrapolation of nucleon masses from the lattice and calculation of the effect 
of the lattice size on nucleon masses are discussed in Section lOl 

If the pions are light compared to the inverse lattice extent, but the lattice size is not 
too small such that 1/L ~ e, where e is a small parameter, one can study fluctuations 
around the zero modes within the so-called e expansion of X-P^ [33] . 



4 LATTICE RESULTS 

4.1 Setting the scale 

Masses and decay constants coming out of a simulation are at first dimensionless numbers 
in units of the lattice spacing. The value of the lattice spacing is determined by calculat- 
ing a suitable quantity aM on the lattice and adjusting the corresponding dimensionful 
quantity M to its physical value. If the calculation is free of systematic errors such as lat- 
tice spacing, finite volume and quenching effects, using any quantity should give the same 
result. In practical calculations all of those errors can occur. Then, 'suitable' means that 
systematic errors of the quantity used to set the scale and the quantity that is supposed to 
be calculated cancel as well as possible. For example, the lattice scale can be determined 
using the static quark potential, which has small discretization errors (O(a^) or higher). 

A typical length scale is tq related to the interquark force [34] with rg ^|^_^ = 1-65, 

which can be calculated on the lattice with high precision [35]. The physical values cor- 
responding to potential models are around tq = 0.49 — 0.5 fm. Unless noted otherwise we 
use ro = 0.5 fm. For the string tension a, usually experimental values of \fa = 427 or 440 
MeV are assumed. Other quantities frequently used to set the scale are the p meson mass 
rUp, the nucleon mass, the decay constants f^^ and fx, and charmonium and bottomonium 
level splittings. 
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Figure 2: Discrepancy of lattice spacings from various physical quantities on quenched 
lattices. Results are from [35,37,39,41,42,44,45]. Average over spin orientations is 
denoted by an overbar. 



np=0 results 











r-e- 


fK 




-e- 


3Mh - Mn 
2Mg - 


-e- 






— e— 




T(1D-18) 




-0- 


T(2P-1S) 




-e- 


T(38-18) 


e 




T(IP-IS) 


1% errs 


1: , , , 





O.B 1.0 1.2 1.4 

LatticeQCD/Exp't 



nF=3 results (2003) 



, , , ;| 






-Bi 




fK 




3M5 - 






2Mb_ - Mt 






i/'(lP-18) 




»- 


T(1D-18) 






T(2P-ig) 






T(38-ia) 


k 




T(lP-lS) 


1% errs i| 




, , 1 , , , , 



0.8 1.0 1.2 1.4 

LatticeQOD/Exp't 



Figure 3: Comparison of lattice with experimental results from [40], using zero (left) and 
Nf = 2iight + Istrange flavors (right), a from the T' — T mass splitting. 



In Fig. 121 we show examples for the discrepancy between lattice spacings from different 
physical quantities in the quenched approximation. 

With two flavors (TV/ = 2), the agreement is improved: using Wilson gauge fields and 
two flavors of 0(a) improved clover sea quarks, Ref. [36] quotes an agreement of scales 
from nip, fx and = 0.5 fm. However, in the two flavor calculations of [37, 38] (Wilson 
gauge fields, staggered sea and clover valence quarks) and of [39,41] (RG gauge fields and 
tadpole-improved clover sea and valence quarks) at a ~ 0.5 GeV~^, a ~ 20% discrepancy 
between lattice spacings from Xb — T mass splittings and rrip remains. 

With two flavors of light and one flavor of strange dynamical quarks, using a 1-loop 
Symanzik O(a^) improved gauge action and a tree- level tadpole 0{a ) improved staggered 
sea quark action Ref. [40] finds an agreement of a variety of physical quantities with 
experiment (see Fig. EI). 
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Figure 4: Unquenched light baryon spectrum from the lattice. Left and middle: results 
from [41]. Right figure: from [42]. Open symbols denote quenched, filled symbols 
unquenched results. K input: strange quark mass set by fixing the K meson mass to 
the physical value, input: strange quark mass set by fixing the (f) meson mass. 



4.2 The light baryon spectrum from the lattice 

In quenched calculations it was found that the features of the experimental light hadron 
spectrum are described well by the lattice [46,47]. It is of interest to study whether un- 
quenching improves the agreement. In Fig.|4]we plot the baryon spectrum from the recent 
unquenched simulations of [41,42]. Discrepancies with experiment of ~ 2a remain. A rea- 
son may be uncertainty in the chiral extrapolation. Ref. [42] assigns additional systematic 
errors of up to 25 MeV from the chiral extrapolation uncertainty and the determination 
of tq. In Table ^ we give results for light baryon mass splittings corresponding to the 
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Table 1: Light baryon mass splittings. The first error is statistical, the second is the 
difference from fixing the strange quark mass using the K or (p meson where applicable. 
The quantity used to fix the lattice scale is indicated. 

results shown in Fig. E) For the splittings, the agreement with experiment is at the l — 2a 
level. A recent calculation [43] with A'^; = 2iight + Istrange finds a A — A^ splitting which 
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agrees well with experiment. 

4.3 Nucleon mass: chiral extrapolation and finite size effects 




Figure 5: Chiral extrapolation of nucleon masses on large lattices in two-flavor QCD using 
non-relativistic (left: from Ref. [48]) and relativistic (right: from Ref. [50]) xPT at 0{p^). 
The dot-dashed line on the left shows the non-relativistic 0{p^) result. 

HBxPT predicts at 0{p^) a correction ~ to the quadratic dependence on ttItt, 
but with a coefficient which is very different from the value found from fits to the lattice 
data. In Ref. [48], a good description of lattice data up to pion masses ~ 600 MeV could 
be achieved using the non-relativistic formalism at 0{p^) (see Fig. on the left). With 
relativistic xPT at O(p^) [49], the agreement with the lattice data is also good up to 
rather large pion masses, as shown in Fig. El on the right [50]. 

Having ensured that relativistic x^T 0{p^) indeed describes the nucleon mass on very 
large lattices, it should be possible to calculate the finite size effects on lattices which are 
not too small within in this formalism. Calculating the difference of the nucleon self- 
energy in a spatially finite and infinite volume within x^T at Olp"^) [50], assuming an 
infinite temporal extent of the lattice, one finds a good agreement with the finite size 
behavior of the lattice results. An example for a pion mass around 550 MeV is given in 
Fig. El 

The non-relativistic formalism at 0{p^) predicts finite size effects which are clearly 
smaller than the finite size effects of the lattice data [51]. 

4.4 The spectrum of hadrons with a b quark 

A heavy quark with infinite mass can be regarded as a color source which is static in the 
rest frame of the hadron and whose spin is not relevant to the interactions. Corrections 
due to the finiteness of the heavy quark mass can be included in a 1/M expansion. Within 
the Heavy Quark Effective Theory (HQET), the mass of a heavy-light hadron H can be 
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Figure 6: Volume dependence of the nucleon mass relativistic X-P^ compared with A'^^ = 2 
lattice results, from Ref. [50]. Solid line: O(p^) result, dashed line: only 0{p^) terms. 



thought of as consisting of the following contributions: 

\h\QD^Q\H) , {H\Qa-BQ\H) 



Mh = Mq + A 



2Mr 



2Mh 2Mh 



0{1/M'), (17) 



where Q is the heavy quark spinor, Mq the heavy quark mass, A the binding energy of the 
meson for Mq — > oo, and the other two terms the expectation values of the heavy quark 
kinetic energy and the spin-colormagnetic interaction energy respectively. We give a brief 
summary of lattice results on the hyperfine splittings B* — B and Bg — Bg in Section l4.4.1[ 
Results on P wave states and baryons are given in Tables |21 El and 0] Masses averaged 
over spin-orientations (spin-averaged) are denoted by an overbar. The first error on the 
individual lattice results includes statistical errors and uncertainties fixing the masses to 
the physical values, the second, where applicable, is a chiral extrapolation uncertainty. To 
calculate weighted averages, we include an estimate of systematic errors from the actions. 
Most of the calculations use NRQCD, except for [52] who uses heavy clover quarks and [53] 
who simulates Bg mesons in the static approximation and interpolates between the static 
and experimental Dg mesons. 

Discretization errors with non-perturbatively 0(a) improved clover light quarks (finer 
lattice of [44] and [55]) are Ola^Aq^j^), whereas the tadpole- improved light clover action 
has O^o^Aq^^) and 0{agaAQCD) errors (coarser lattice of [44] and [56,57]). Refs. [58,59] 
use O(a^) tree- level tadpole-improved clover light actions respectively. Ref. [60] uses 
staggered light quarks. The scale has been set with rrip except in Ref. [57] sets the scale 
with ^/a = 427 MeV. Ref. [53] and [54] uses tq with physical values of 0.5 and 0.525 
fm respectively, and [60] uses tq = 0.5 fm and quarkonia at and around the charm [61]. 
All other calculations from the set discussed here use rrip. The NRQCD action has errors 
O^asAqcD/ M) from corrections to the spin-magnetic coefficient. Errors on spin splittings 
are treated as being dominated by an error on the spin- magnetic coefficient of 0{ag) ~ 
20 - 30%. 
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The systematical error of each result is divided into a part common to all calcula- 
tions, which is taken to be of the order of the error of the calculation with the smallest 
uncertainty, and and a rest which is treated as independent. The error on the average is 
rescaled by r = a/xVC^^'^^^^' where N is the number of results, if r > 1. The second 
error on the averages comes from the variation due to the 10% ambiguity between using a 
from nip and a from tq = 0.5 fm in the quenched case, and asymmetric chiral extrapola- 
tions where applicable. The Xb mass difference is not included in the estimate of the 
scale variation since it gives values for spin-independent mass splittings which are much 
higher than experiment. For example, Ref. [38] quotes Bg — B and the Afy — B splittings 
of 118 and 670 MeV if the scale is set with the Xb — splitting instead with rrip. The 
experimental values are 90 and 345 MeV respectively [62]. 

If a collaboration quotes results from several lattice spacings, they are plotted starting 
from the coarsest lattice on the left. Asymmetric errors are added linearly in the plots. 

For the error estimates we use nominal values of Aqcd = 400 MeV, M = 5 GeV and 
as = av{l/a), where is defined in the potential scheme described in [63] at the 

scale q*. Since the lattice results have rather varying central values we do not calculate 
the error in percent of the individual lattice splittings but of the experimental splittings 
or nominal estimates thereof. 

4.4.1 B* -B splitting 

Results from quenched lattice NRQCD calculations of the B* — B and B* — Bg splittings 
fixing the scale with nip (e.g. [44,56,57,59,64]), and tq = 0.5 fm [60] are found to be 
around 25 — 35 MeV, compared to experimental values of 45.8(4) MeV and 47.0(2.6) 
MeV [62], respectively. Using relativistic 0(a) improved heavy quarks, Ref. [45] obtains 
a splitting around 10 — 20 MeV with an error of 10 MeV. A quenched calculation using 
the FNAL action [65] setting a with the charmonium IP — IS splitting quotes a. B* — Bg 
splitting of around 40(10) MeV. A preliminary calculation comparing results with zero 
and two flavors of 3x the strange quark mass [57] on lattices with a ~ 0.2 fm finds an 
unquenched value of around 33 MeV with only an insignificant increase compared to the 
quenched result from coarse lattices. A recent NRQCD calculation with 2-1-1 dynamical 
flavors [66] finds a B* - Bg splitting of 42.5(3.7) MeV. 

4.4.2 Orbitally excited B mesons 

For heavy-light mesons it seems appropriate to use a hydrogen-like picture for the coupling 
of angular momenta of the quarks. In the infinite mass limit, there are two P wave energy 
levels with light quark angular momentum ji = 1/2 and 3/2. At finite M there is an 
additional hyperfine structure due to the coupling of the heavy quark spin. This results 
in one level with angular momentum zero (Bq), two with angular momentum one {B[ and 
Bi), and one with angular momentum two {B2). Here we discuss the Bq and i?2 level 
splittings. 

Lattice results are presented in Table |21 In the error estimation, we use guessed values 
where no experimental value is available: 400 MeV for Bq— B and B*q — Bg and 500 MeV 
for B^- B and 5*2 - Bg. 

Experimental knowledge of the P state level structure is still sparse. The Particle 
Data Book lists two candidates. There is the 5} (5732) resonance with a width of ~ 130 
MeV, which is believed to come from several narrow and broad P wave states, and the 
-B*j (5850) signal with a width of ~ 50 MeV, which can be interpreted as stemming from 
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Figure 7: Comparison of splittings between P wave and the ground state of Bg mesons 
from Refs. [44,53,56,58,60] (lattice) and Refs. [67,69,70] (models). The lines denote the 
experimental value of the narrow (5850) resonance believed to come from orbitally 
excited Bg mesons [62]. 



excited B^ states. A comparison of the lattice results with experiment and with model 
calculations is given in Fig. [7| 

The sign of the -83 —-Bq niass difference is disputed among potential model calculations 
(e.g. [67-70]). Individual lattice calculations [44,58] find a splitting around zero and are 
within errors compatible with a small negative splitting, but the lattice average for i?2 — -Bq 
is positive. 

4.4.3 b baryons 

Baryons with one b quark can be thought of as two light quarks coupling to form a spin 
zero or spin one diquark. The state with a spin zero diquark is the A;,. If the diquark has 
spin one, the heavy quark can couple to a spin 1/2 state, the T^b, and a spin 3/2 state, 
the Eft*. If the light quarks in the T,b and the are substituted by strange quarks one 
obtains the Qb and the Qb- 

In Table m we summarize results for the spin-independent splittings Ab — B and E5 — A^ 
with M(Sb) = [2M(Sfe) +4M(Sb*)]/6 and M(B) = [3M(S*) + M(5)]/4, and the E§ - Sf, 
hyperfine splitting. The values quoted for [55] are obtained by interpolating the heavy- 
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Figure 8: — B splitting from quenched (left) and Nj = 2 (middle) lattices (Refs. [38, 
52,55-57,59]). The lattice average shown on the right is quenched. CP-Q and CP-2 
denote quenched and unquenched results from [57], respectively. Model results are from 
Refs. [71,72]. 



strange meson mass to the Bs mass and interpolating the baryon splittings to the thus 
obtained b quark mass. 

We compare lattice results with calculations within a constituent quark model [71] and 
a Skyrme model [72] where the baryon is described as a bound state between a soliton 
and a heavy quark. 

The experimental values of the Ab — B and the Ac — D splittings are very close: 311(9) 
and 310(2) MeV respectively. The only 1/M correction to these splittings comes from the 
heavy quark kinetic energy. Its contribution to the A;, mass appears to be very close to 
the contribution to the B mass. 

The quenched lattice average differs from the experimental value by less than 2a. 
Preliminary results with two flavors of dynamical quarks around the strange quark 
mass [38, 57] are even higher. Ref. [57] finds an increase of ~ 15% if Nf is changed 
from zero to two dynamical quarks of around 3x the strange quark mass. Calculation at 
smaller quark masses would clarify whether the difference is related to a chiral extrapo- 
lation uncertainty. 

In Fig. El we show results for the spin- averaged Ef, — Af, splitting. Ref. [72] gives a rela- 
tion between the heavy-light and light baryon splittings, AM (Sq — Aq) /AM (A — A^) = 

2/3, for Q = c,b. For Q = c, the equality holds well experimentally. Lattice results for 
the ratio with Q = b vary between 0.5 and 1. 

In Fig.iniwe also give results for the Eft — E^ hyperfine splitting. The expectation from 
HQET is that the hyperfine splitting is generated by the spin-chromomagnetic interaction 
(see Eq. (fTTjl ) and should be proportional to 1/Mq. Rescaling the experimental value of 
AM(E* — Ec) — 67 MeV by the ratio of b and c quark masses, one expects AM(Efe — E^) ~ 
20 MeV, and the lattice results are compatible with this expectation. The error on the 
Eft — Eft and fib — ^6 splittings is estimated to be 0{as) x 20 MeV. 
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Figure 9: Above: S;, — A;, splitting from the lattice (left) and model calculations (right). 
The dashed line shows the experimental value for the Sc — A^ splitting. Below: S^* — S;, 
splitting. Lattice data from Refs. [52,55,56,59], model results from Refs. [71,72]. 



4.5 fB 

In Table El we summarize lattice results for fs and /^^ since 1998. The first error given 
in the Table is statistical, and the second is the systematical error given by the authors 
added in quadrature. 

First we address unquenching effects on Jb- They depend on how the scale is set. In 
Table El we compare ratios of decay constants from quenched and two-flavor simulations 
with the same gauge field and valence quark actions and find an increase of ~ 10% if a is 
set with /tt, 10 — 20% if a is set with ntp and no increase with T and (for fo^) with tq. 

We calculate weighted averages for quenched, Nf = 2 and Nf = 2 + 1 results of fs 
and fs^- Since the methods for error estimation can vary considerably between different 
collaborations, even if similar lattice actions and parameters are used, we make new 
assignments motivated by the error analysis of the authors themselves. We assign common 
systematic errors to the calculations with RG gauge fields using NRQCD [73] and using 
the FNAL heavy quark action as non-relativistic effective field theory without taking 
the continuum limit [74], and to the calculations using Wilson gauge fields and NRQCD 
(quenched are [64,75,77]). For quenched configurations with clover light quarks, we assign 
systematic errors of 20 and 22 MeV respectively for Jb and Jb,- Ref. [60] uses NRQCD 
with an O(a^) tadpole improved gauge action and staggered light valence quarks, and we 
use their own systematic error assignment. 

The quenched calculations of Refs. [45, 78-80] use Wilson gauge and clover quarks 
at a = 0.35 — 0.37 GeV~^ (/? = 6.2) simulated at the charm quark mass. Ref. [79] 
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uses tree- level tadpole-improved clover quarks without including 0{as x a) terms in the 
renormalization. Ref. [78] uses a non-perturbatively 0(a) improved clover quark action 
and a partly non-perturbative current renormalization. Refs. [45,80] use nonperturbative 
0{a) improvement except for a perturbative value for the O(asamq) quark mass correction 
to the renormalization constant. Although different degrees of improvement are used, and 
the scaling behaviour is found to be different, the results for fr,^ agree at (3 = 6.2. We 
therefore assign a common systematic error to these results. According to the estimate 
of the discretization error given in [79] (8%) and of a 1/M extrapolation error of ~ 9% 
given in [45] we use 23 MeV for and 26 MeV for f Be- 
liefs. [76,81,82] use heavy quarks in the FNAL formalism and extrapolate their results 
to a — >■ 0. Refs. [83,84] use a step scaling method with the Schrodinger functional and 
clover heavy quarks. Part of the renormalization factors is calculated nonperturbatively. 
Their results are continuum extrapolated. Ref. [85] uses an interpolation between static 
and clover charm quarks which are non-perturbatively improved using the Schrodinger 
functional and continuum extrapolated. 

The second error on the quenched results includes the ambiguity between scales from 
rrip and T level splittings by varying the result by -|-30% if the scale is taken from rup or 
^/a = 427 MeV, (-3 + 27)% if the scale is set with and (-12 + 18)% if the scale is set 
with ro = 0.5 fm. The lattice spacings calculated in [79, 80] from K physics are close to 
the results using /,r at the same /3 values, and the scales from Jk determined in [81] are 
close to the ones using nip from the same actions and P values. We also quote the upper 
bound from the variation between scales from rrip and tq = 0.5 fm (second error given 
in square brackets) which is slightly lower than but within errors compatible with the 
unquenched central value with Nf = 2 + 1. The uncertainty in fixing the strange quark 
mass is also included. Where it is not given by the authors, we include a (+7) MeV error. 

For two-flavor QCD, we estimate the systematic errors from non-relativistic meth- 
ods [36,38,73,74] to be 21 MeV for and 25 MeV for /b^. For the systematic error 
of the continuum extrapolated results with FNAL heavy quarks [82,88] we use the esti- 
mate of Ref. [82] including errors due to continuum extrapolation, perturbation theory, 
1/M extrapolation and, where applicable, the spin-magnetic coefficient. Ref. [36] makes 
an estimate of the light quark mass dependence of fs using 1-loop xPT and quotes an 
uncertainty (^g) MeV on Jb from the chiral extrapolation. We assign the same error also 
to the Nf = 2 results of Refs. [38,73,74,88]. For the values of Ref. [82] we use their own 
estimate of the chiral extrapolation error. The uncertainty in the chiral extrapolations, 
a ~ 50 MeV increase in the decay constants if the Xb — splitting instead of rup is used 
to set the scale quoted by [38] and [73], and the variation from determining the strange 
quark mass by setting the K or (p meson mass to the physical value give the second error 
on the Nf = 2 averages. The variation between using rrip and to set the scale and 
the chiral extrapolation uncertainty determined by [82], and the variation in the strange 
quark mass give the second error in square brackets. 

The results of this procedure are included in Table El as averagel. 
Refs. [74] and [38] quote statistical errors on their quenched and Nf = 2 results on 
respectively which are about half the statistical errors of other calculations with 
a similar ensemble size. We therefore also calculate the average with their statistical 
errors and the error of Ref. [84] enlarged by a factor of two. The result, which is very 
close to averagel, is quoted as average2 in Table El Except for the uncertainties due to 
fixing the scale and reaching the physical quark masses, the errors on these averages are 
only few percent. The double ratio of B and D decay constants /b^ ^J MbJ {fsy/MB x 
foy/M]:)/ (fn^ ^JMd^) should be independent of the exact form of the chiral extrapolation 
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up to 1/Mq corrections, as is supported by the results of a two-flavor calculation using 
FNAL heavy quarks [88]. As argued by Ref. [87] using x^T, the chiral extrapolation 
uncertainty of the ratio /^^ ^J MbJ {fB\/ Mb) x fn/ fx should also be small. 

Employing unquenched staggered {Nf = 2+1) MILC gauge field ensembles at a ~ 0.13 
fm, the NRQCD estimate of [86] is: 

/b. = 260(7)(28), (18) 

where systematical errors are added in quadrature. Calculations of the decay constants 
with the staggered iV/ = 2 + 1 configurations using NRQCD [89] and FNAL [90] heavy 
quarks are in further progress. Within the statistical and systematical errors quoted in 
Tabled the results with A'^^ = 0, 2 and 2 + 1 agree among each other. 

We relate this to the experimental value for ffj^ using unquenched lattice results for 
the ratio fBs/fos from two-flavor calculations which work directly at the b and c quark 
masses without using extrapolations. 

Taking the experimental value fo^ = 283(45) MeV (Eq. ©), and the range of values 
for the ratio fBslfos from Table El one obtains /^^ = 230 — 260 MeV. 

Other recent review articles [92-96] quote lattice estimates for Jb and which are 
within errors in agreement with the averages quoted in Table 

In Table El we compare the lattice results with recent sum rule [97-100] and potential 
model calculations [101], and we find that they are within errors in agreement. 



5 Conclusions 

Applications of non-relativistic QCD and chiral perturbation theory in lattice calculations 
are presented. The status of lattice results on the light and heavy-light hadron spectrum 
and the decay constants and /^^ is summarized, and weighted lattice averages for h 
hadron mass splittings and decay constants are calculated. The agreement of the hadron 
spectrum with experiment is a major success of lattice QCD in general, and of non- 
relativistic methods for heavy quarks in particular, and supports the reliability of lattice 
predictions of hadronic matrix elements. The lattice has become instrumental in QCD 
calculations. Work on further understanding and reduction of lattice errors is in progress 
and will enable very precise checks. 
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Ref. AM{B) MeVJ AM{Bs) MeV 


B*-B 


Lattice 


[44 
[44 
[58 
[56 
[60 
[60 
[60 
[60 
[60 
[60 
[60 
[53 
[54 
ave 


, a ~ l.lGeV^^ 
, a ~ 2.6GeV-i 

, 11 
, 11 
, 11 
, 11 

, Asq 
, Asq 
, Asq 

rage 


400(30)(-) 
295(45)(ii) 

mil) mil) 

374(37) 357(27) (°) 
442(56) 
471(25) 
315(105) 
425(60) 
285(78) 
523(94) 
403(56) 
386(31) 

408(67) 419(37) 
402(41)(f) 382(23)(?D 


Model calculations 


[67 
[69 
[70 




450 440 
453 466 
427 431 


B*-B 


Lattice 


[44 
[44 
[58 
[56 
[57 
[53 
[54 
ave 


, a ~ l.lGeV-^ 
, a ~ 2.6GeV"i 

rage 


402(78)(-) 
474(62) (J6) 

493(1) 478(i) 
526(45) 493(26) 
426(17) 

534(52) 

440(77) 455(41) 
498(49)(f) 487(31)(r) 


Model calculations 


[67J 
[69] 
[70] 


450 490 
453 469 

435 447 


Preliminary experimental 


62 




419(8) 483(15) 



Table 2: B orbital excitations. Only statistical errors and, where quoted by the authors, 
errors due to chiral extrapolation and fixing the b quark mass are shown. For the B2 
of [44] at a ~ l.lGeV"^ we quote the result from the lattice operator they use for the 
P — S splitting. The other lattice operator corresponding to the B2 in their calculation 
gives a ~ 45 MeV higher result. 
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Ref. 


AM{B) 


MeV 


AM{Bs) 


MeV 












* 










Lattice 




44 
44 
58 
56 
53 
54 
ave 


, a ~ l.lGeV"^ 
, a ~ 2.6GeV~i 

rage 


18(1) 
155(32) 

32(87) 
98{47)0 


41(94)(-) 
179(65)(^) 

27(i) 

136(23) 

148(61) 

36(55) 

101(25)(ii) 




Model calculations 


[67 
[68 
[69 
[70 




40 

-155 

-5 
7 




50 

3 

16 







Table 3: P state fine structure of B mesons. Only statistical errors and errors due to 
cliiral extrapolation and fixing the b quark mass are shown. 
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Ref. Ab-B[MeV] - A^MeV] 


Lattice 




52 
56 
59 
59 
55 
57 
38 
57 
que 


, as = l.lGeV-i 
, as = 0.9GeV-i 

, quenched 
,Nf = 2 
,Nf = 2 
nched average 


338(-) 
370(67) 

mz) 

367(}?g) 
389(44) 
361(22) 
545(40) (22) 
417(19) 
372(33) (f) 


186(-) 

221(71) 

156(1) 

191(1?) 
122(65) 

174(28) (f) 


Models 


[72J 
[71] 

[71]c quark 


312 


196 
217 
212 


Jl/xperiment 


62 


311(10) 




Ref. - Sfe [Me V ill - ih MeV 


Lattice 


[56J 

[59], as = l.lGeV-i 
[59], as = 0.9GeV-i 
[55] 

average 


19(7) 
22(12) 

nil) 

10(12) 
18(8)(^) 


18(4) 

18(1) 
20(9) 
7(4) 
13(7)(g) 


Models 




71 

72 




10 

8 





Table 4: b baryons. Only statistical errors and, where applicable, systematical errors due 
to scale setting except for the quenched scale ambiguity and fitting are shown. 
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Ref. 




scale 


/B[MeV 


/i?.[MeV] 


Lattice 








N 


/ = o 






76] 
81] 
77] 
64] 




nip 
Ik 

TTl 

= 427MeV 


173(4)(13) 

147(ll)(li) 
170(5)(15) 


199(3)(14) 
185(^3)(9) 
175(8)(16) 
191(4)(17) 




78] 
45] 
79] 
75] 
74] 
"73] 




fn 
h 

nip 
nip 
nip 


173(13)(f) 

195(6)(i) 

177(17)(22) 

188(3)(26) 
191(4)(27) 


196(ll)(o ) 
220(6)(i) 
204(12)(i) 
187(4)(15) 
220(2) (i) 

220(4)(31) 




[80] 
82] 
83] 
84] 
85] 
[60] 

average 1 
average2 


U 
ro 
ro 
ro 
ro 


174(22) (8) 

173(6)(16) 
170(11) (23) 


204(15)(8) 

199(5)(i) 
192(9)(25) 
192(6)(4) 
205(12) 
225(9) (34) 

198(5)(-[^1) 

201(6)(S^''^) 








N 


/ = 2 






38] 
74] 
73] 
82] 
36] 
88] 




nip 
nip 
rrip 

U 
nip 

ro = 0.49fm 


186(5)(25) 
208(10) (29) 
204(8) (29) 
190(7)(??) 
191(10)(^°) 
181(7)(i^) 


215(3)(-) 
250(10)(|) 
242(9) (34) 
217(6)(i) 
215(9)(}|) 


average 1 




190(10) (^^''') 


226(15) (f^^) 


average2 






226(15)(f^^) 










= 2 + 1 






86] 


T(2S' - IS) 


260(7)(28) 


Sum rules 


[97J 
[98] 
[99] 
[100] 




203(23) 
206(20) 
210(19) 
180 - 190(30) 


236(30) 
244(21) 


Potential models 




101 






178(15) 


196(20) 



Table 5: fs and /^^ from the lattice. Statistical errors and systematical errors given by 
the authors are included, adding the systematical errors in quadrature. The method to 
set the scale is indicated in the second column. The first error on the averages is due to 
the statistical and systematical errors of the individual results, while the second error is 
from chiral extrapolation uncertainties and scale ambiguities as explained in the text. 
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Ratios of decay 
Ref. scale fB^~'^/ fB^~^ 


constants 

Ref. fnJfD. 


[82 
[82 
[74 
[73 
[73 
[91 


fn 1.10(6) 

rUp 1.19(6) 
rUp 1.11(6) 
TTip 1.07(5) 
T{P-S) 0.97(5) 
(/z,J ro 0.98(4) 


[74J(iVj = 0) 0.88(1) 
[82](7V^ = 0) 0.891(12)(^0) 
[85] {Nf = 0) 0.81(6) 
[74](A^; = 2) 0.94(6) 
[82](iV; = 2) 0.922(13)(i) 



Table 6: Ratios of decay constants. The first error is statistical, the second the system- 
atical errors given by the authors added in quadrature, where applicable. 
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